Ferromagnetic-like closure domains in ferroelectric ultrathin films. 
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We simulate from first-principles the energetic, structural, and electronic properties of ferroelectric 
domains in ultrathin capacitors made of a few unit cells of BaTiOa between two metallic SrRuC>3 
electrodes in short circuit. The domains are stabilized down to two unit cells, adopting the form of a 
domain of closure, common in ferromagnets but only recently detected experimentally in ferroelectric 
thin films. The domains are closed by the in-plane relaxation of the atoms in the first SrO layer of 
the electrode, that behaves more like SrO in highly polarizable SrTiOa than in metallic SrRuC>3. 
Even if small, these lateral displacements are essential to stabilize the domains, and might provide 
some hints to explain why some systems break into domains while others remain in a monodomain 
configuration. An analysis of the electrostatic potential reveals preferential points of pinning for 
charged defects at the ferroelectric-electrode interface, possibly playing a major role in films fatigue. 

PACS numbers: 77.80.Dj,77.22.Ej,77.84.Dy,68.55.-a 



Ultrathin film ferroelectric capacitors are under active 
investigation [H, Q . Of considerable technological inter- 
est as memories, transducers, and electromechanical de- 
vices, they present problems of considerable scientific in- 
terest. Although technologically relevant films are thicker 
than 100 nm, deeper understanding of the origin of these 
problems requires combined experimental and theoreti- 
cal studies of thinner regimes. On the one hand, recent 
breakthroughs on materials synthesis and characteriza- 
tion techniques have allowed the growth of ferroelectric 
thin films with a control at the atomic scale and the local 
measurement of the ferroelectric properties 0]. On the 
other, the steady increase in computational power and 
improvements in the efficiency of the algorithms permit 
accurate first-principles study of larger and more complex 
systems, overlapping in size with those grown epitaxially. 

Prominent among problems of interest is understand- 
ing the mechanisms screening charge densities at the in- 
terfaces. The termination of the ferroelectric polariza- 
tion at the surface or the electrode interface generates 
a polarization charge which gives rise to a depolarizing 
field tending to suppress the polarization. Two mecha- 
nisms are traditionally invoked for the compensation of 
the polarization charges: the first, screening by charge 
accumulation at the electrode (or even by ionic adsor- 
bates 0, Q); the second, the breaking up of the system 
into domains @, 0] • 

Previous first principles local density calculations on 
realistic short-circuited ferroelectric capacitors suggested 
that a monodomain configuration for the polarization 
was unstable below a critical thickness that ranged be- 
tween m — 2 and m = 6 layers [1, 0, [l(| of ferroelectric, 
depending on the perovskite, the electrode, and the ter- 
mination at the interface. In all these approaches the 
electrode was the only source of screening, providing free 
charges that accumulate at the interface on the metallic 
side and even decay exponentially into the first few layers 
of the ferroelectric, or sharing the ionic displacements re- 



sponsible for the polarization in the ferroelectric ll| . In 
any case, the mechanism is ineffective below this critical 
thickness where the paraelectric phase was stabilized. 

In this letter we simulate from first principles, 
within the local density approximation to the den- 
sity functional theory and the numerical atomic orbital 
method as implemented in the Siesta code 12 1, typi- 
cal SrRuOs/BaTiOs/SrRuOs ferroelectric capacitors, in 
which we allow the system to form domains. Our start- 
ing point is the reference paraelectric heterostructure de- 
scribed in Ref. § , that is now replicated N x times along 
the [100] direction, where N x ranges from 2 to 8. A soft 
mode distortion of the bulk tetragonal phase is super- 
impossed to the BaTiOa layers of the previous paraelec- 
tric configuration, so the polarization points upwards in 
half of the superlattice and downwards in the other half 
[see inset of Fig. QJa)]. The twinning on both the BaO 
(Ba-centred) , and Ti02 (Ti-centred) planes is considered. 
Then, the atomic positions of all the ions, both in the 
electrode and in the ferroelectric thin film, are relaxed 
till the maximum component of the force on any atom 
is smaller than 0.01 eV/A for m = 2, and 0.04 eV/A 
for ttl — 4. Very accurate computations are required 
since the differences in energy between relevant phases 
are eight orders of magnitude smaller than the absolute 
value of the energy. The electronic density, Hartree, and 
exchange-correlation potentials are computed in a uni- 
form real space grid, with an equivalent plane- wave cut- 
off of 400 Ry. Once self-consistency is achieved, the grid 
is refined (reducing the distance between points by half) 
to compute the total energy and atomic forces. We used 
a Afc x x 12 x 1 Monkhorst-Pack mesh for all the Brillouin 
zone integrations, where = 4p- except for the inter- 
face with N x —8, where Nk x = 2. All the calculations 
are performed at T = 0. Details on pseudopotentials and 
basis set used can be found in Ref. fliij ]. 

Our calculations support stabilization of a polydomain 
phase with an exceptionally small periodicity below the 
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previous critical thickness [see Fig. QJa)], in good a gree - 
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ment with the results obtained with Landau theory [14 1 . 
For a two unit cell thick film, m = 2, the extra source 
of screening is efficient provided that the domain period 
is between two and four times the thickness of the film. 
Within this region, the energy cost of forming the do- 
main wall is compensated by reduction of the net po- 
larization charge at the interfaces. As in 180° stripe- 
domains in bulk [l5j ]. the Ba-centred wall configuration 
is preferred. The energy difference between the most sta- 
ble polydomain and the paraelectric phase for a capacitor 
with to = 2 is very small, of the order of 1.5 meV (~ 16 
K) for the whole supercell. For this thickness there is 
essentially no energy difference between domains of lat- 
eral periods N x = 4 and 6, suggesting that both might be 
equally present in a sample. Heating or cooling processes 
might help the system to overcome potential energy bar- 
riers and activate the transition between them. Although 
the conductive nature of the substrate is different, this 
fact might provide an extra source of explanation [16| for 
the intriguing richness in behaviour of the stripe domain 
patterns observed experimentally in PbTiC>3 thin films 
rown on SrTiC>3, where two different periods coexisted 
(Note that our ratio between domain periods, 1.5, 
is close to the experimental factor 1.4 for the so-called a 
and phases in Ref. 0.) 

The energy differences between polydomain and para- 
electric phases increase very quickly with thickness [Fig. 
Gib)] and amounts to 120 (80) meV for a Ba-centred (Ti- 
centred) domain wall capacitor with N x = 4 and to = 4. 
For this size, the polydomain phases are more stable than 
the monodomain configuration, itself more stable than 
the paraelectric phase by 20 meV. 

The minimum energy structures of these ferroelectric 
capacitors, shown in Fig. [2] display the closure domain 
configuration proposed by Landau and Lifshitz 17| and 
Kittel for magnetic systems. At the centre of the 
BaTiC>3 layer, the displacement of the atoms and there- 
fore the corresponding local dipolcs, point normal to the 
interface (coordinate z), as expected for 180° stripe do- 
mains. However, approaching the ferroelectric/electrode 
interface a small tilt towards [100] is observed. Remark- 
ably, the domains are not closed by the surface layer of 
the ferroelectric [3] but by the in-plane displacements 
of the Sr and O atoms at the first layer of the electrode, 
which yield a closure domain pattern, with 90° domain 
walls with the z oriented domains inside the film. In 
contrast to the metallic relaxations in monodomain con- 
figurations, where ionic displacements penetrate into the 
metal over a distance of two or three unit cells 
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the displacements beyond the second R11O2 layer are 
negligible, an indication of more effective screening pro- 
duced by the domains of closure. The in-plane displace- 
ments of the atoms at the interfacial SrO layer, although 
small in magnitude, stabilize the domain structure. If 
a constrained relaxation is performed in which the in- 



0.4 
0.3 
0.2 





SrR 


U0 3 




—I 


Bal 


10, J 






SrR 








SrTiOj 





m=2 






D 








M -A 

IN — *f 






X 




n 
u 










-5 


A 

A 


> 






E, 


-10 


■ ' 






• 


z* 
>, 


-15 










Enei 


-20 


□ ' 




-25 






-30 









-35 





FIG. 1: (color online). Difference in energy between polydo- 
main and paraelectric phases as a function of (a) the domain 
period N x for a ferroelectric thin film two unit cells thick 
(m = 2), and (b) the thickness of the ferroelectric film for 
a capacitor with N x — 4. The energy of the paraelectric 
phase (dotted line) is taken as reference. First-principles re- 
sults for both Ba-centred (circles, solid line) and Ti-centred 
(squares, dashed) domain walls are shown. In (a) differences 
in energies between local minima of the polydomain phase 
are represented by error bars. Inset, structure of the ferro- 
electric capacitor considered. N x is the stripe period and m is 
the thickness of the ferroelectric thin film, in number of unit 
cells of the ferroelectric perovskite oxide. In (b) the result for 
the most stable monodomain configuration is also shown (tri- 
angle). Full symbols correspond to constrained relaxations 
where no in-plane displacements are allowed. 



plane forces on all the atoms are artificially eliminated, 
the atoms move back to the paraelectric positions for 
to = 2, or to a structure comparable in energy to the 
most stable monodomain configuration for m = 4 [Fig. 
[ijh)]. Whether the in-plane displacement is allowed or 
not might partially explain for the very different configu- 
rations found experimentally in related heterostructures: 
Lichtensteiger et al., using the same experimental setup, 
have observed how high-quality ultrathin films of PbTiC>3 
grown on Nb-SrTiC>3 electrodes remain in a monodomain 
configuration [2l| (although with reduced polarization 
and tetragonality) whereas they form domains when the 
electrode is replaced by Lao. 67Sro.33Mn03 1241 ■ The same 
domain formation is suggested for Pb(Zro.2Tio.s)03 on 
SrRu0 3 [13. 

Regarding the origin of this polarization induced relax- 
ation, the analysis of the projected density of states (not 
shown here) shows that the SrO layer closest to the inter- 
face behaves more like SrO in SrTi03 than SrO in metal- 
lic SrRu03. Similar behaviour was found in AO/ATi03 
heterostructures 13], where A = Ba or Ti. Both first- 
principles computations [24| and experimental measure- 
ments [25J have shown that SrTiOs is highly polarizable 
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when combined with BaTiOs in heterostructures. 

Similar domain patterns have been found using a first- 
principles effective hamiltonian for Pb(Zr .4Ti .6)O3 [f6| 
asymmetrically screened (grown on a nonconducting sub- 
strate and with a metal with a dead layer as top elec- 
trode), and using a Landau- Ginzburg phenomenological 
approach for a PbTi03 thin film [261 ]. both asymmetri- 
cally and symmetrically coated with insulating SrTi03. 
Here, the domains of closure are obtained even for a sym- 
metrical metal/fcrroelectric/metal capacitor, where the 
metallic plates should provide significant screening. 
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FIG. 2: (color online). Schematic representation of the 
atomic relaxations in patterns of domains of closure with do- 
main period of N x =4 (a), and N x = 6 (b). Balls, represent- 
ing atoms, are located at the positions of the reference para- 
electric phase. Atomic displacements for the polydomain con- 
figuration after relaxation are represented by arrows, whose 
magnitude can be gauged with respect to the displacements in 
the bulk tetragonal phase of BaTiOa at the scale on the left. 
Dotted lines indicate the position of the domain wall. Only 
Ba-centred domains are shown. Similar results are obtained 
for Ti-centred domains. 



The polarization can be estimated from the structural 
calculations. Figure [3] displays how much the polar dis- 
tortion along z is changed by the presence of a domain 
pattern. We define as A the average of the change of 
distance, with respect the most stable paraelectric con- 
figuration, between a Ti atom and the nearest O atom 



lying on top along the z direction (cf. Ref. 15), nor- 
malized with respect to the short Ti-0 distance in the 
tetragonal bulk phase. A norm is a very sensitive indica- 
tor of the polar order: it is zero as long as the atoms lie 
in the paraelectric position and tends to unity as the full 
bulk polar distortion is attained. Figure [3] shows a very 
narrow 180° domain wall, about a lattice constant wide, 
across which the polar distortion symmetrically reverses 
its sign. In contrast to 180° domains in bulk [15j, where 
the ferroelectric distortion fully recovers its bulk value 



by the second atomic plane far away from the domain 
wall, here it only amounts to 13% of the bulk value at 
the centre of each domain for a m — 2 structure, suggest- 
ing that the polarization for the thin film is one order of 
magnitude smaller than in bulk. This mean polarization 
increases with thickness, and already amounts to 60% for 
a thin film four unit cells thick (m = 4). 
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FIG. 3: Measurement of the polarization in the ferroelectric 
layer as a function of position along the [100] direction of the 
capacitor, (a) Definition of the average change in distance A 
between Ti and apical O in a chain along [001] for an interface 
with m = 2. In every case, the atomic positions correspond 
to the lowest energy structure. A positive value of A means 
a polarization pointing upward, (b) Profile of the normalized 
averaged change in distance along z as a function of the posi- 
tion of the chain for a Ba-centred interface of domain period 
N x = 4. The chains are numbered as indicated in Fig. [2] 
Results are shown for m = 2 (dashed line) and m = 4 (dot- 
dashed). Dotted lines represent the position of the domain 
walls. 



Ideally, closure domains do not produce any polariza- 
tion charge since the normal component of the polariza- 
tion is preserved across any domain wall. Therefore the 
depolarizing field should vanish everywhere (l8j . and a 
constant electrostatic potential is expected. To further 
check this point we plot in Fig. [5] the nanosmoothed 
27, [28| electrostatic potential along z as a function of 
the position along the [100] direction of the capacitor. 
No nanosmoothing is performed along x. For a stripe 
of thickness m = 2 and period = 4, the potential is 
esentially flat at the centre of the domain, in contrast to 
the depolarizing field reported for monodomain configu- 
rations [Ij]. A large microscopic field along [100] appears 
inside the domains of closure at the metal-ferroelectric 
interface. The origin of this field is due to the difference 
in polarization in the domain of closure (the polarization 
along x equals 0.9 and 3.3 /zC/cm 2 for m = 2 and m = 
4, respectively) and inside the thin film in our realistic 
capacitor. Besides, after nanosmoothing in z a residual 
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depolarizing field along [001] is identified in the neigh- 
borhood of the domain wall, decaying rapidly away from 
it. This last field might be responsible for the lowering 
of the polarization with respect to bulk shown in Fig. [3] 
Both fields might play an important role in the fatigue of 
ferroelectric capacitors, the most serious device problem 
in ferroelectric thin films [lj ■ In particular we identify at 
the ferroelectric/electrode interface the preferred points 
of migration of charged defects, which pin the domain 
walls and inhibit their motion [29[ . The depolarizing field 
at the centre of the domain increases with the domain pe- 
riod; it starts to be appreciable for N x = 6 [Fig. HJb)] 
and finally desestabilizes the ferroelectric distortions for 
N x = 8, as shown in Fig. QJa). 
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FIG. 4: (color online). Map of the nanosmoothed electro- 
static potential in a two unit cells thick ferroelectric capaci- 
tor with a stripe period of N x =4 (a), and N x — 6 (b). The 
arrows represent the atomic displacements with respect the 
paraelectric phase as in Fig. [2] Only the displacements of the 
cations are shown for simplicity. Full lines are a schematic 
representation of the domains of closure, while dashed lines 
mark the position of the BaTiOs/SrRuOs interface. 



Although we have demonstrated that the domains are 
stable, it is not clear whether the capacitor as a whole 
can be called ferroelectric since, for this, the polarization 
has to be switchable under external electric fields [ 1 9l. 1 2 31] - 

Our calculations provide insightful results on the 
energetic, structural and electronic properties of 
ferromagnetic-like closure domains in ultrathin capac- 
itors, only recently observed experimentally by J. F. 
Scott's group [13]. We provide some hints to explain why 
some systems break into domains while others remain in 
a monodomain configuration. We also predict the pref- 
erential sites for pinning charged defects, important for 
understanding the fatigue of thin films. 

We thank M. H. Cohen for the critical reading of the 
manuscript. Calculations were performed on the com- 
puters of the ATC group and on the Altamira Supercom- 
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